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Groundstates of certain materials can support exotic excitations with a charge 
that's a fraction of the fundamental electron charge. The condensation of these 
fractionalized particles has been predicted to drive novel quantum phase tran- 
sitions, which haven't yet been observed in realistic systems. Through numeri- 
cal and theoretical analysis of a physical model of interacting lattice bosons, we 
establish the existence of such an exotic critical point, called XY*. We measure 
a highly non-classical critical exponent rj = 1.49(2), and construct a universal 
scaling function of winding number distributions that directly demonstrates 
the distinct topological sectors of an emergent Z 2 gauge field. The universal 
quantities used to establish this exotic transition can be used to detect other 
fractionalized quantum critical points in future model and material systems. 
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It is a remarkable fact that, in this age of high-energy accelerator experiments, certain types 
of fundamental quantum particles can only be studied in what seems the relatively pedestrian 
world of tabletop condensed-matter physics experiments. Consider for example the familiar 
electron, carrying fundamental charge e. Unlike the proton, whose charge originates from 
quarks with fractional charge, no energy is sufficiently high to break up the charge of an elec- 
tron. However, as demonstrated years ago by the measurement of fractional Hall conductance, 
if one places an electron in certain clean two-dimensional materials in a strong magnetic field, 
its charge can indeed break into fractions - e/3, e/5, and so on (1) - each fractional charge 
arising from a quasiparticle emerging in the sample. 

Such quasiparticles share all the important characteristics of real particles. In "deconfined" 
phases of matter with a gap to excitations, the quasiparticles can be separated a large distance 
from each other, making them well-defined localized objects with a sharp energy-momentum 
dispersion relation. Many types of quasiparticles are known to occur as elementary excitations 
in the low-temperature groundstates of quantum materials, most typically as the more mundane 
examples of normal bosons or fermions. u The condition for their emergence is the realization 
of a particular exotic "vacuum", which is the groundstate of some condensed matter system. It 
is as excitations out of this strange vacuum that these unfamiliar quasiparticles can flit in and 
out of existence. 

For over a decade, condensed matter physicists have searched for delicate and elusive frac- 
tionalized particles in systems other than Hall effect materials. Theoretical predictions have 
identified a class of low-temperature paramagnets, the quantum spin liquids, as holding partic- 
ular promise for supporting them (2-4). However, experimental searches for these fractional 
charges and their parent spin liquid vacuum remain unconvincing. This is due in part to the 
difficulty in constructing measurements that are able to identify the experimental signatures 
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necessary to betray their existence in the variety of specific materials - ranging from fabricated 
solid state devices, to delicate organic magnets - thought to harbor possible spin liquid states (5). 
Perhaps the most interesting recent experimental candidate is a set of materials that may dis- 
play fractional particles with a gapless Fermi surface (6-9). This leads to an even more tenuous 
situation: fractionalized quasiparticles without an excitation gap, where interaction between 
quasiparticles makes it problematic even to define a fractionalized excitation. □ In bosonic 
systems, such interacting gapless fractional quasiparticles have been proposed to mediate novel 
quantum critical points that exist in certain order-to-order transitions (10). These "Landau- 
violating" critical points (which rely on long-wavelength fluctuations of fractional particles) 
have the advantage that signatures of the fractionalization are manifest in universal quantities, 
such as critical exponents, avoiding reliance on measurements of specific material-dependent 
quantities. 

A major goal of the theoretical community has been to demonstrate the existence of these 
"deconfined" quantum critical points in realistic microscopic models, a task that necessarily 
falls to large-scale numerical simulation - quantum Monte Carlo (QMC). A variety of models 
which may contain such order-to-order critical points (the best candidate being Sandvik's J-Q 
model (11)) are currently under intense numerical study. In this paper, we study deconfined 
quantum criticality in a different context - an order-to-disorder transition between a superfluid 
and a gapped spin liquid state, in a physical model of lattice bosons. We establish the existence 
of fractionalization at this critical point by confirming the prediction (12-15) of a strongly 
non-classical universal critical exponent rj. As suggested in Ref. (10), the defining property of 
interacting fractionalized excitations is the presence of an emergent gauge field. By considering 
topological properties of this gauge field, which enforce nonlocal constraints, we construct 
the first direct measurement of fractionalized excitations at a critical point. These topological 



2 However, in the case of Refs. (6-9), Fermi liquid theory suggests that fractionalized excitations remain sharp 
near the surface. 



Figure 1: A section of the toroidal kagome-lattice simulation cell. Gold spheres label sites 
occupied by bosons in a representative configuration. Red spheres are fractional charges, mark- 
ing "defect hexagons" (those which do not have three bosons per site). Defect hexagons are 
shown joined by an open string. We also show a representative closed string not associated with 
fractionalized particles. 

properties manifest themselves in new universal scaling functions that can be used as a test of 
the existence of fractionalized excitations in a variety of other models and materials. 

The Kagome Bose-Hubbard Model. We examine a variant of the Bose-Hubbard Hamil- 
tonian, introduced in Ref. (16), which is a simple model of hard-core bosons hopping (with 
strength t) on a two-dimensional kagome lattice with an energetic constraint (V) favoring three 
bosons per hexagonal lattice plaquette: H = —tJ2{ij) [K^j + hb]] + V J2o( n o) 2 - This model 
has been convincingly shown to contain a superfluid to spin-liquid quantum phase transition at 
7/i ~ 7 through a smoking-gun measurement of the topological entanglement entropy (16). 
In that work, two universal exponents were measured at the critical point, related to the diver- 
gence of the correlation length (u) and the isotropy of the space and time dimensions (z), both 
of which fall in the conventional 3D XY universality class. A section of the toroidal kagome 
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lattice geometry used in our QMC simulations, along with a representative boson configuration, 
is shown in Fig. [TJ 

This system conserves total charge (the boson particle number). Deep in the spin liquid 
phase (V 3> t), it is known that exotic quasi-particles carrying half-odd integer charge exist 
as fundamental excitations out of the groundstate vacuum. At infinite V and at half filling, the 
system is restricted to a space of states with three particles in each hexagon. In this limit the 
model can be mapped to a dimer model on a triangular lattice formed by the centers of the 
hexagons, with three dimers per site (77). Turning on a non- vanishing t/V produces virtual 
excitations out of this space of states, leading to effective exchange terms which give a quantum 
spin liquid phase. Imagine removing a single particle from this state, creating a pair of defect 
sites in the dimer model with only two dimers. In a phase with mobile, deconfined defects, these 
defects can become separated by large distances, and since the state has total charge —1 relative 
to the vacuum, each defect carries charge —1/2. Similarly, defects with four dimers carry 
charge +1/2. The model is thus in a Z 2 spin liquid phase, but electric defects carry half-odd 
charge^ Overlaying this configuration with a fixed reference configuration gives a transition 
graph, which is a gas of strings, with open strings connecting defects (18). There are three 
dimers per site so strings may intersect; in Fig. [T] we only show two strings for clarity. Since 
links in the string arise alternately from the given configuration and the reference configuration, 
links need not have a particle on them in the image shown. 

Colloquially, one can view the fractionalized particles as a "square-root" of the boson field, 
with the Z 2 gauge field present due to the sign ambiguity in the square-root. Let , b denote 
creation and annihilation operators for the real physical bosons, and let <f>\ (ft denote the same 
operators for the fractionalized particles. At an XY* critical point, the fractionalized field 



3 This can be viewed as a topological quantum field theory, with particles 1, e +1 / 2 , e^ 1 / 2 , m, where e ±x / 2 
denotes the charge on the electric particle and the ... denote an infinite sequence of additional particles with charge 
shifted by any integer. 
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Figure 2: Equal time Green's function G(L/Q) = (tf(0, 0)fe(L/6, 0)) as a function of system 
size L. 

4> undergoes an ordinary XY transition. As pointed out in Refs. (14, 19), this dramatically 
affects the observed critical exponents, since the order parameter b usually associated with the 
transition is actually a composite operator, made up of two (f> fields. This leads to the same 
exponent v as in the ordinary XY critical point. However, the exponent rj controlling the equal 
time ground state correlation function (tf(0)b(x)) is significantly modified, leading to 77 ~ 1.45 
as estimated by field theoretic and Monte Carlo simulations of the correlation of a composite 
operator in the 3D XY model (20, 21). This contrasts strikingly with 77 fa 0.03 in the ordinary 
3D XY transition. 

The measurement of 77 in the kagome model is challenging, since QMC does not directly 
access the superfluid order parameter operators necessary to measure the relevant correlation 
function. To overcome this fact, we have measured the equal time Green's function in real 
space: G(x) = G(t = 0,r) = (tf(0)b(r)} (22,23). This measurement involves keeping track 
of the defects created in the non-local loop algorithm (24) as it traverses the QMC's d+ 1 space- 
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time simulation cell. At the critical point, G(r) should decay at sufficiently long distances as 
l/r 1+TJ . Finite size effects can be minimized by looking at G(L/m), which decays as 1/L 1+V , 
where L is the linear system size and m is some (fixed) number. In Fig. [2l we show G(L/6) 
as a function of L. The Green's function decays algebraically with 77 = 1.49(2), a value that is 
consistent with 77 for a composite operator in the 3D XY model. This strongly non-classical 77 
shows that the boson operator b is indeed a composite operator of fractionalized excitations. 

More direct evidence can be obtained from topological properties, as we now show by con- 
structing new scaling functions that exploit universality, and as such can be used as a test for 
the existence of fractionalized excitations in other models and experiments. Consider a path in 
imaginary time in which a total of n x of the <\> particles wind around the x-direction of the torus 
geometry. In this case, the total charge W x that winds around the torus is equal only to n x /2. 
However, the only closed paths in imaginary time that are allowed are those in which n x is even, 
as the physical charge W x winding around the torus must be an integer. This is a manifestation 
of topological properties of the Z 2 gauge field. The gauge field is gapped, so does not affect the 
dynamics of the </> field, but there are a total of eight different topological sectors (two sectors in 
each space direction and one in the time direction) and summing over sectors restricts to even 
winding and even number of the fractionalized particles (see Fig. [Q where winding a particle 
drags a string, changing the topological sector). It is known that the winding number variance 
(W%) in the ordinary XY model is a universal function of L/v/3, where L is the length scale of 
the torus and v is a non-universal velocity. It is natural (and we have confirmed numerically) that 
the full winding number distribution, as quantified by the probabilities P(W X ) of observing a 
given winding, is also a universal function. Thus, the universal critical properties of the winding 
number distribution W x in our model can be computed from a model undergoing an ordinary 
XY transition in 2 + 1 dimensions, by restricting to sectors with even winding numbers n x , n y 
and even total particle number. Then, one can determine the probability of observing given n x 
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Figure 3: Comparison of winding number distribution P{W X ) for the kagome-lattice XY* 
transition at L = 16, 24, to P even (2W x ) for a conventional L = 16 square-lattice transition 
(described in the Supporting Online Material, Section 1.1). Dashed lines show the unprojected 
distribution of XY winding numbers P(W) for comparison. Upper curves are W = 1 and 
lower are W = 2. The x-axis is fi/L for the kagome lattice, and (3/L times a non-universal 
velocity ratio for the XY model. 



where P even (2W x ) is the probability, after projecting onto these even sectors, of observing n x = 
2W X . While one does not know the scaling functions controlling P even analytically, we have 
sampled them using quantum Monte Carlo on a different model which contains a conventional 
XY transition (Supporting Online Material, Section 1.1). As shown in Fig. [3] the prediction of 
the theory shows excellent agreement for two distinct scaling functions (we also tested P(3) but 
it is not shown) over the majority of the observed curve. The region in which the functions do 
not agree, at higher temperatures, is due to low energy magnetic excitations (Supporting Online 
Material, Section 1.2); one can observe that as system size is increased, the range of agreement 



with 



P(W X ) = P even (2W x ), 



(1) 
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This provides a direct test of the fractionalized charge. The fit of all the scaling functions 
uses only a single adjustable parameter, the non-universal velocity ratio vxy /vkagome- In the 
limit ofv/3/L — > oo, the winding numbers become large, and both the projected and unprojected 
winding number distributions converge to approximate Gaussians, so the projection onto even 
numbers has little effect on the variance. Then, one will observe {W%)xy* = 0-/ty{W%)xY, 
where the subscripts denote the two different universality classes and we assume both mod- 
els are at the same ratio v/3/L. However, testing this agreement is difficult as it requires an 
accurate measurement of the non-universal velocity. In fact, the best way we have found to 
measure the non-universal velocity (as the high temperature measurements are complicated by 
the magnetic excitations) is to study the full distribution of scaling functions which provides a 
much more stringent test. Note also that the scaling functions of the unprojected XY model are 
very different from the projected functions at moderate v/L/3, where both functions are strongly 
non-Gaussian and non-Poissonian. 

The Phase Diagram. As mentioned above, Fig. [3] shows slow convergence to universal 
scaling in the high temperature critical regime (small v/3/L). To understand this, we have 
considered the phase diagram of the model as a function of V/t and temperature T. At small 
V/t, the model is superfluid at zero temperature, with a second order transition at non-zero 
temperature. At larger V/t, but still below the critical point, the phase transition becomes first 
order, similar to behavior seen in Ref. (25). This is likely due to vortices where the phase of 
winds by %, which are more relevant at long length scales than the usual 2n vortices. Close 
to the critical point, the transition should return to being second order at non-zero temperature. 
However, the presence of low energy Z 2 magnetic vortices has a strong effect on the non-zero 
temperature properties. 

The scenario that is consistent with the numerical data is that for V ^> t, the electric defects 
have an energy of order V and so are of much higher energy than the magnetic defects. Increas- 
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ing t reduces the energy of the electric defects due to gain in kinetic energy, until this energy 
vanishes at the transition. However, even a few magnetic defects can strongly influence the 
winding number distribution; the number of magnetic defects is roughly L 2 exp(— const. /T), 
and so for T ~ the number of magnetic defects is exponentially suppressed, but for mod- 
erate L, the prefactor leads to a substantial effect manifest as slow convergence of the scaling 
function at small 0/L as shaded in orange in Fig. [3] 

Discussion. We have conclusively established the existence of a phase transition with de- 
confined fractionalized excitations in quantum Monte Carlo simulations of a physical model 
of lattice Bosons. The topological properties of this transition are demonstrated by universal 
scaling functions of the winding number. This test of the fractionalized charge uses the full 
winding distribution, and may be regarded as a strongly interacting analogue of current fluctu- 
ation measurements used experimentally to test fractionalized charge, as in Ref. (26). The full 
winding distribution is in principle experimentally accessible, as it is the Fourier transform of 
the free energy as a function of flux through the torus. 

One can imagine more general deconfined quantum critical points with other gapped gauge 
fields, such as a Zk gauge field. Lattice models for such theories are lacking, but if found, the 
winding number distribution that we have introduced will present a clear test of the universality 
class. The fact that topological properties remain important even at the critical point suggests 
that such systems, despite being gapless, might be useful for quantum computing applications 
in the future. 

References and Notes 

1. R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

2. D. S. Rokhsar, S. A. Kivelson, Phys. Rev. Lett. 61, 2376 (1988). 



10 



3. M. Hermele, et al, Phys. Rev. B 70, 214437 (2004). 

4. Y. Tang, A. W. Sandvik, arXiv: 1107.1439 (201 1). 

5. L. Balents, Nature 464, 199 (2010). 

6. M. Yamashita, et al, Science 328, 1246 (2010). 

7. M. S. Block, D. N. Sheng, O. I. Motrunich, M. P. A. Fisher, Phys. Rev. Lett. 106, 157202 
(2011). 

8. Y. Shimizu, et al., Phys. Rev. Lett. 91, 107001 (2003). 

9. Y. Kurosaki, et al., Phys. Rev. Lett. 95, 177001 (2005). 

10. T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, M. P. A. Fisher, Science 303, 1490 
(2004). 

11. A. W. Sandvik, Phys. Rev. Lett. 98, 227202 (2007). 

12. A. V. Chubukov, T. Senthil, S. Sachdev, Phys. Rev. Lett. 72, 2089 (1994). 

13. A. V. Chubukov, S. Sachdev, T. Senthil, Nuclear Physics B 426, 601 (1994). 

14. T. Graver, T. Senthil, Phys. Rev. B 81, 205102 (2010). 

15. S. V. Isakov, T. Senthil, Y. B. Kim, Phys. Rev. B 72, 174417 (2005). 

16. S. V. Isakov, M. B. Hastings, R. G. Melko, Nature Physics, doi:10.1038/nphys2036 (201 1). 

17. L. Balents, M. P. A. Fisher, S. M. Girvin, Phys. Rev. B 65, 224412 (2002). 

18. D. S. Rokhsar, S. A. Kivelson, Phys. Rev. Lett. 61, 2376 (1988). 

19. T. Senthil, O. Motrunich, Phys. Rev. B 66, 205104 (2002). 

11 



20. P. Calabrese, A. Pelissetto, E. Vicari, Phys. Rev. E 65, 0461 15 (2002). 

21. H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, A. M. Sudupe, Phys. Lett B 387, 125 
(1996). 

22. N. V. Prokof'ev, B. V. Svistunov, I. S. Tupitsyn, Physics Letters A 238, 253 (1998). 

23. A. Dorneich, M. Troyer, Phys. Rev. E 64, 066701 (2001). 

24. O. F. Syljuasen, A. W. Sandvik, Phys. Rev. E 66, 046701 (2002). 

25. S. V. Isakov, A. Paramekanti, Y. B. Kim, Phys. Rev. B 76, 224431 (2007). 

26. A. Bid, N. Ofek, M. Heiblum, V. Umansky, D. Mahalu, Phys. Rev. Lett. 103, 236802 (2009). 

27. C. Castelnovo, C. Chamon, Phys. Rev. B 76, 184442 (2007). 

28. S. Iblisdir, D. Perez-Garcia, M. Aguado, J. Pachos, Phys. Rev. B 79, 134303 (2009). 

29. R. K. Kaul, R. Melko, Phys. Rev. B 78, 014417 (2008). 

Acknowledgments The authors thank T. Senthil for enlightening discussions. MBH thanks 
the Aspen Center for Physics for hospitality. This work has been supported by the Natural 
Sciences and Engineering Research Council of Canada (NSERC) and the Swiss HP 2 C initiative. 
Simulations were performed on the Brutus cluster at ETH Zurich and the computing facilities 
of SHARCNET. 



12 



1 Supporting Online Material 
1.1 XY Model Simulations 

In this section, we describe quantum Monte Carlo (QMC) simulations of the two-dimensional 
XY model with an applied staggered field: 



where Si is a spin- 1/2 operator, and X\ and y-i are the lattice coordinates of each spin. This 
Hamiltonian was chosen as a simple quantum lattice model to realize a transition in the XY 
universality class between superfluid and insulating phases. It has XY symmetry but no higher 
continuous symmetry (in contrast, the XX Z model has SU (2) symmetry at the Heisenberg 
point), and applying a staggered field rather than a uniform field was chosen to avoid producing 
a chemical potential for bosonic excitations of the model which would produce terms linear in 
the time derivative for the effective action of the bosonic field. We employ the stochastic series 
expansion QMC method with directed loops (the code was checked on a 4 x 4 lattice against 
Lanczos diagonalization), to study the phase transition out of the superfluid phase with increas- 
ing h s /J. The phase transition is identified through finite-size scaling of the superfluid density, 
p s , as illustrated in Figure |4j The broad features of the superfluid density identify a transition 
near h s ~ 1 (inset), which through a finite-size scaling collapse is seen to be consistent with 
a continuous quantum phase transition with exponents z = 1 and v = 0.6717 - the values for 
the 3D XY universality class. With the systems sizes studied, the best collapse of the finite-size 
data occurs at a critical h s /J = 0.995(5). 

Simulations were run to extract the probability distribution of the winding number squared 
at h s /J = 0.995. Measurements of both the full distribution were considered, as well as the 
"projected" distribution. In the latter case, data was recorded only when both (x and y) winding 
number were even (or zero), and the particle number was even. These distributions allow us to 




(2) 



13 



1.2 
1 

Pi^0.8 
0.6 
0.4 
0.2 

-4 -2 2 4 

[(V/)-(V-/)^ 1/v 

Figure 4: Finite size scaling collapse of the superfluid density of the XY model in a staggered 
magnetic field, using v = 0.6717 and (h s /J) c = 0.995. 

calculate both the superfluid density of the physical bosons, as well as the "projected" space of 
composite bosons, which simulates the composite bosons of the fractionalized model (Fig. [5]). 

We have tested the universality of the winding number distribution, to verify the hypothe- 
sis that at the critical point this converges to a universal function of v/3/L. In Fig. 0, we plot a 
scaling collapse of the projected winding distribution functions P even (21V) as a function of j3/L 
for L = 16, 20, 24. Only very minor differences can be seen between the collapsed curves; as 
a further test we plot the same winding distributions but allow additional non-universal factors 
(of the order of 2%) to multiply the x-axis for each system size. As can be seen in Fig. [5l the 
collapse can be made almost perfect with these factors, indicating that the dominant finite size 
effect is a slight L-dependent correction to the velocity. Since our comparison to the kagome 
lattice model involves multiplying by a non-universal velocity ratio, such L-dependent correc- 
tions to the velocity are unimportant in comparing the projected XY and kagome curves. 

1.2 Low Energy Magnetic Excitations 

The theory of the XY* critical point predicts that the electrical particles, carrying ±1/2 charge, 
become gapless as one approaches the critical point from the Z 2 spin liquid. Indeed, in order 
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Figure 5: Left: Plot of projected winding numbers, P even (2W) as a function of (3/ L for L = 
16, 20, 24. Right: Plot of projected winding numbers, P even (2W) as a function of /3 for L = 
16, 20, 24, with small non-universal scaling applied to a>axis for each size. Upper curves are 
W = 1 and lower curves are W = 2. 

to have the correlation function (&t(0)fe(r)) not decay exponentially, some excitation carrying 
charge must become gapless. The theory also predicts that magnetic particles remain gapped 
as the critical point is approached. Once the electric particles condense in the superfluid phase, 
the magnetic excitations become confined and vanish from the spectrum of the theory. 

At a non-zero temperature, however, there is some density of thermally excited magnetic 
particles in the spin liquid phase and at the critical point. In this subsection we consider the 
properties of these excitations. In particular, one complication is that while their excitation gap 
is predicted to remain non-zero at the transition, the actual numerical value is quite small, mean- 
ing that they can have a significant effect even at what seems to be a relatively low temperature. 

First, consider the Z 2 spin liquid. The study of (16) showed two plateaus in the topological 
entanglement entropy as a function of temperature. These plateaus are characteristic of two 
different energy scales, one for electric and one for magnetic excitations (27, 28). Deep in the 
Z 2 spin liquid side with f < V, clearly the lower energy scale corresponds to magnetic exci- 
tations, that energy scale being set by t 2 /V, while the electric excitation energy scale is set by 
V . However, as t/V is varied, the evidence suggests that the characteristic energy scale of the 
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Figure 6: Schematic phase diagram of the kagome Bose-Hubbard model. 



magnetic excitations remains quite small. This implies that at fairly small temperature there still 
is a non-zero density of magnetic excitations. However, close to the critical point, we are con- 
fronted with an apparent paradox: the upper plateau appears to be at a fairly high temperature, 
and appears to correspond to the electric excitations. One piece of evidence for this is based on 
the fact that the fluctuations in particle number are small (and correspondingly the uniform sus- 
ceptibility described below is small) at this temperature scale, and also the fact that the density 
of defect hexagons becomes small at this scale. However, the energy for the electric excitations 
must go to zero as they are the excitations that carry non-zero charge. The resolution of this 
paradox that we propose consistent with the data is that the energy of the electric excitations 
does go to zero, but that at a small non-zero temperature there is also a density of magnetic ex- 
citations. These magnetic excitations renormalize the energy of the electric excitations upwards 
(the energy of an electric excitation can be thought of as a balance between a positive term due 
to creating a defect hexagon and a negative term due to gain in kinetic energy, but the magnetic 
excitations reduce the kinetic energy gain), and so we see a regime at intermediate temperature 
which has a large density of magnetic excitations but a low density of electric excitations. This 
regime is similar to the cooperative paramagnet phase (25); it is not separated from the high 
temperature paramagnet by a phase transition but there is a very distinct, measurable crossover 



16 



in properties such as the topological entanglement entropy and the susceptibility. 

These magnetic excitations also have an important effect on the properties for V/t < (V/t) c , 
in the superfluid phase. Consider the phase diagram as a function of V/t and temperature, Fig.|6] 
At small V/t we see a continuous phase transition as a function of temperature. This transition 
is likely Kosterlitz-Thouless, but closer to the critical point the transition appears to become first 
order (see also Ref. (25)). As discussed in the text, this is likely due to the existence of vortices 
where the phase of the field winds by ix. However, such phase winding corresponds to binding 
a Z-2 vortex to the phase vortex and hence the importance of these vortices depends upon the 
energy scale of the magnetic excitations. Since the T = phase transition is continuous as a 
function of V/t, eventually we expect that for V/t sufficiently close to (V/t) c the transition as 
a function of T will become continuous again, but we have not fully resolved this region. 

Evidence for the existence of the first order transition can be obtained, for example, by 
studying the uniform susceptibility as a function of V/t and temperature. This susceptibility is 
obtained from fluctuations in the particle number by 

Xu = j^((n-n) 2 ), (3) 

where n is the total boson occupation number. As shown in Fig. [71 we see a rise in susceptibility 
at low temperature, followed by a sudden decline. We associate this with a first order transition 
to a cooperative paramagnet. More direct evidence is a visible double peaked structure of the 
distribution of the kinetic energy, E k = —t(J2ij(.b]bj + fejfrj)) - see Fig.[8l 

The existence of these magnetic excitations poses one further issue for analyzing the critical 
point. One might hope to follow a similar approach to Ref. (29), where the uniform suscepti- 
bility and superfluid density are both considered as function of /3/L, and universal amplitudes 
are extracted from considering the combination of these two quantities. The superfluid density 
becomes non-negligible once (3/L is large, while the uniform susceptibility is non-negligible 
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Figure 7: The uniform susceptibility Xu as a function of temperature for different values of V/ 1 
on the superfluid side of transition. Upper panel: L = 8. Middle panel: L = 16. Lower panel: 
L = 24. 
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Figure 8: Distribution of the kinetic energy per site close to the finite-temperature transition at 

V/t = 6.5 and (3 = 7.5. 

when /3/L is small (since the gap to a particle excitation is of order 1/L, the particle number 
variance is exponentially small in (3/L). If the limit of large L is taken at fixed (3/L, then we 
also take the large (3 limit and hence we will eventually reach a temperature scale low enough 
to have a negligible number of magnetic excitations. Note that the total number of magnetic 
excitations is roughly L 2 exp(— (3AE), where AE is the energy gap to magnetic excitations, 
and hence this number decays exponentially in L assuming fixed (3/L. However, in the regime 
where the uniform susceptibility is non-negligible, which corresponds to a small value of (3/L, 
we find that there still is a non-negligible number of magnetic excitations, due to the L 2 pref- 
actor, at accessible values of L. For this reason, as shown in Fig. |9l we do not see a good 
scaling collapse at small and intermediate (3/L. One may expect that a scaling collapse might 
be recovered in some range at large values of (3/L for large enough system sizes L, when the 
gap to electric excitations becomes well bellow the gap to magnetic excitations. Some evidence 
is shown in Fig. [71 where one can see that the uniform susceptibility is strongly suppressed for 
small system sizes but it becomes finite for large enough system sizes as the two energy scales 
cross over. The effect gets more pronounced as we approach the quantum critical point. 
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Figure 9: Superfiuid density p s and the uniform susceptibility Xu as a function of LT. The 
uniform susceptibility is strongly suppressed at low and intermediate temperatures. Lines guide 
the eye. 

While this means that we cannot find universal properties by analyzing x M and p s in par- 
allel, the winding number distribution measurement discussed in the main text provides a way 
of extracting the universal properties of the fractional charge by studying properties that are 
determined by properties at a sufficiently low temperature to remove the magnetic excitations. 
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